The main purpose of this article is to increase the efficiency of the least squares method in numerical solution of Ill-posed functional and physical equations. Determining the least squares of a given function in an arbitrary set is often an Ill-posed problem. In this article, by defining artificial constraint and using Lagrange multipliers method, the attempt is to turn n-dimensional least squares problems into (n-1) ones, in a way that the condition number of the corresponding system with(n-1) -dimensional problem wIll be low. At first, the new method is introduced for2 and 3-term basis, then the presented method is generalized for n-term basis. Finally, the numerical solution of some Ill-posed problems like Fredholm integral equations of the first kind and singularly perturbed linear Fredholm integral equations of the second kind are approximated by chain least squares method. Numerical comparisons indicate that the chain least squares method yields accurate and stable approximations in many cases.